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I. INTRODUCTION 


Nowadays, we are witnessing an outstanding progress in observational cosmology, espe¬ 
cially owing to the precise measurements of the cosmic microwave background (CMIU, which 
is one of the best known observational windows to the early epochs of the universe [1| . Such 
a breakthrough provides the opportunity for theoretical physicists to test the predictions 
of their theories about those early stages. In particular, the possibility that we may falsify 
our hypotheses about the origin of the universe in this way, makes clear the necessity of 
providing a quantum theory for the gravitational interaction, as General Relativity suffers 
a predictability breakdown in this regime: The big-bang singularity. 

One of the most promising approaches to accomplish such a quantization of gravity is 
based on the formalism of loop quantum gravity (LQG) Q. The application of its techniques 
to describe the quantization of cosmological models leads to the branch of research known as 
loop quantum cosmology (LQG) [^. LQG has proven to produce some remarkable results in 
the analysis of homogeneous models. In particular, the initial big-bang singularity is avoided 
and replaced with a quantum mechanism, called the big-bounce [^, |^. However, testing the 
robustness of LQG calls for obtaining physical results from more realistic scenarios, such as 
inhomogeneous cosmologies. In this context, a hybrid approach has been proposed for the 
quantization of models of this kind, based on the assumption that the most relevant quantum 
geometry effects would mainly affect the homogeneous degrees of freedom. Such assump¬ 
tion involves a splitting of the phase space of the system into a homogeneous sector and an 
inhomogeneous one. Then, a loop quantization is adopted for the homogeneous degrees of 
freedom (which, in this way, fully retain the genuine quantum features of the space-time), 
while the inhomogeneous degrees of freedom are treated by means of a more conventional 
Fock quantization. This quantization strategy was applied for the hrst time to the case of 
the Gowdy cosmologies with linear polarization of the gravitational waves and with the spa¬ 
tial topology of a three-torus, T^, achieving a complete quantization of the model jsl, [^. The 
study of these tractable cosmologies provides the opportunity to develop approximate meth¬ 
ods and techniques to solve the complicated dynamics of inhomogeneous systems. These 
methods could become particularly interesting when it comes to analyzing the dynamical 
behavior of more realistic models such as Friedmann-Robertson-Walkeiij (FRW) geometries 
with cosmological perturbations, in the context of an inflationary universe. Recently, these 
systems have drawn a substantial attention within the framework of LQG. They are being 
analyzed by means of this hybrid approach [^, as well as employing other related, but dif¬ 
ferent strategies (see jsl for a quantum treatment of inhomogeneities over a “dressed metric” 
that accounts for a quantum background without any back-reaction, or for an effective 
point of view arising from the requirement that the algebra of constraints closes). 

In order to pursue further the analysis of inhomogeneous cosmologies in LQG, we keep 
investigating the dynamics of the hybrid quantum Gowdy model with linear polarization 
and T^-topology. Specihcally we will consider the system with local rotational symmetry 
(LRS), consistent in this model because it possesses two axial Killing vectors which are 
indistinguishable in principle, and only one direction is anisotropic (namely, the direction 
of the inhomogeneities, in which the gravitational waves vary). Besides, to include matter 
in the system, we will consider a minimally coupled massless scalar held with the same 


symmetries as the metric [1^. The classical phase space of this Gowdy model can then be 


seen as that of a hat anisotropic LRS Bianchi I model with inhomogeneities (gravitational 


^ These cosmologies are also called Friedmann-Lemaitre-Robertson-Walker cosmologies by many authors. 





3 


waves and matter field) propagating in the anisotropic direction. After a partial gauge 
fixing, the reduced system is still subject to two global constraints: The zero-mode of the 
Hamiltonian constraint and a momentum constraint for the inhomogeneous helds. The 
Hamiltonian constraint operator resulting from the hybrid quantization of this model has a 
rather complicated action on the states of the kinematical Hilbert space. Thus, approximate 
methods were developed in ll| in order to find families of states that were approximate 
solutions to this Hamiltonian constraint. These solutions were actually shown to satisfy 
as well, within certain approximations, a Hamiltonian constraint that could be seen as 
corresponding to an FRW universe coupled to a homogeneous massless scalar field and with 
a perfect fluid. In this paper we will further generalize this family of states so that, while still 
being approximate solutions to the full Gowdy model, in addition they approximately obey 
the dynamics corresponding to an FRW cosmology with the kind of geometrical corrections 
that one would expect to find in modified theories of gravity 


The Hamiltonian constraint operator of our hybrid Gowdy cosmology involves two trou¬ 
blesome contributions obstructing its resolution. The first of them is an anisotropy term 
that contains what can be seen as the momentum of the Bianchi I anisotropy variable. It 
acts as a difference operator, coupling the isotropic part of the homogeneous sector with 
the anisotropies. The second contribution is a term that couples the interaction between 
the inhomogeneous modes (which come from both the gravitational waves and the massless 
scalar held) with the homogeneous sector. The discussion in [Tl| and the generalization 
presented here provide quantum states on which the action of these two complicated terms 
can be approximated and disregarded when compared to the other terms in the Hamiltonian 
constraint. Specihcally, in the homogeneous sector (which is characterized by the volume of 
the Bianchi I background and the variable accounting for the anisotropy), the dependence of 
these states on the anisotropy is given by some particular Gaussian-like prohles. Essentially, 
these profiles are sharply peaked at a large value of the anisotropy variable, while being 
reasonably centered at a vanishing value of its momentum. The analyses in [ll| dealt with 
the cases where this peak is either a constant or a volume-dependent function, showing that 
these states provide approximate solutions to the Gowdy model that dynamically behave 
as those of an FRW cosmology coupled to a homogeneous massless scalar field (when the 
peak is a constant) and possibly with a perfect fluid (if the peak depends appropriately on 
the volume), including the case corresponding to a cosmological constant. The extension 
of this family of states that we address in this paper considers the more general case in 
which the peak may depend on a rather generic operator of the homogeneous and isotropic 
geometry, not given necessarily by a function of the homogeneous volume. It will be argued 
that, by imposing certain restrictions on this operator, these states provide again approxi¬ 
mate solutions to the full Gowdy model and, in turn, to the Hamiltonian constraint of an 
FRW universe coupled to a homogeneous massless scalar held, filled with a perfect fluid (as 
in 0 ). and now also geometrically corrected by homogeneous and isotropic curvature-like 
terms. As anticipated above, such new types of corrections can be seen as those expected 
to arise from certain modified theories of gravity, such as f{R) theories (l^ . 


It is worth noting that the considered family of quantum states is not peaked at all at 
homogeneous and isotropic trajectories, but rather the opposite. Nevertheless, in the dy¬ 
namics they behave as if they were homogeneous and isotropic, as much as the Hamiltonian 
constraint is concerned, a fact that can be understood as a quantum collective behavior of 
their anisotropies and inhomogeneities. Let us also comment that this family is not as lim¬ 
ited as one might have thought in principle: Just on the contrary, there exists a considerable 
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variety in the geometry operators that may be used in the construction of the mentioned 
Gaussian-like prohles. 

The organization of the paper is as follows. In Sec. [TTl we will revisit the results of 
the hybrid quantization of the Gowdy model with T^-topology, linear polarization, LRS, 
and minimally coupled to a massless scalar held [l^. In Sec. UTIlwe will discuss how, by 
considering the mentioned Gaussian-like states (already introduced in 0 , and generalized 
here), the complicated Hamiltonian constraint operator can be approximated by a much 
simpler one when acting on them. In Sec. IIVI we will justify under which conditions these 
states can be seen as approximate solutions of the Gowdy model that behave dynamically as 
those corresponding to a modihed FRW cosmology, along the lines commented above. In Sec. 
IVl we will conclude by summarizing our results and discussing their possible implications. 
Finally, we have included two appendices with extra details of the calculations. 


II. HYBRID QUANTIZATION OF THE GOWDY MODEL 

Let us consider the Gowdy model with linear polarization, LRS, and minimally coupled 
to a massless scalar held 4) with the same symmetries as the metric. We will denote the three 
orthogonal spatial coordinates as 6, a, and 6, each of them dehned on the circle. This choice 
of coordinates is adapted to the symmetries of the system, so that the two inhomogeneous 
helds (matter and gravitational waves) have spatial dependence only in e.g. 6. We can then 
expand these helds in Fourier modes in this coordinate. The reduced phase space resulting 
from a partial gauge hxing can be split in two sectors: A homogeneous one formed by 
the zero-modes of the Fourier expansion, which can be identihed with the phase space of a 
LRS Bianchi I system minimally coupled to a homogeneous massless scalar held 0 (the zero¬ 
mode of $); and an inhomogeneous sector containing the non-zero Fourier modes of both 
the linearly polarized gravitational waves and the matter held, as well as their canonically 
conjugate momenta. The whole of this reduced system is subject to two global constraints: 
A momentum constraint, Ce, that generates rigid rotations in 6, and the zero-mode of the 
Hamiltonian constraint, Cq, that generates time reparameterizations. This last constraint 
consists of two terms, namely Cq = Gbi + Ginh- Here, Gbi coincides with the Hamiltonian 
constraint of the LRS Bianchi I model coupled to the homogeneous massless scalar 0, and 
Ginh rules the dynamics of the inhomogeneities, coupling them with the homogeneous sector. 

We then follow the hybrid approach for the quantization of this model [1^. We adopt a 



A. Loop quantization of the homogeneous sector 


In order to adopt a loop representation of the LRS Bianchi I phase space (which is four 
dimensional), it is convenient to introduce four specihc classical real variables to describe 
it, related with the non-vanishing components of the SU{2) Ashtekar-Barbero connection 
and with the non-vanishing components of the densitized triad (see, e.g., 0). Following 


the conventions and notations of |10l. Ill[|. we denote these variables as v, b, A^, and be. The 


absolute value of v is proportional to the volume of the Bianchi I universe (homogeneous 
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sector of the Gowdy model) and Xg measures the anisotropy in the ^-direction. From their 
Poisson brackets (where h denotes the reduced Planck constant), 


{A0,n} = O, {6,n} = ^, {bg,Xg}=‘^^, 

{Ae,6} = 0, {bg,v}=‘^, {bg,b} = ^{bg-b), 


( 2 . 1 ) 


one observes that this set of variables is not canonical. However, it is a convenient choice 
to formulate the improved dynamics scheme in the loop representation 1^. The Hilbert 


space resulting from this quantization is the completion of the linear span of eigenstates of 
the operators v and Xg with respect to a discrete inner product. Explicitly, an orthonormal 
basis of the homogeneous gravitational sector is {|n,A6)) = |n) 0 lA^)} (with v,Xg G M), 
where {v',Xg\v,Xg) = Let us note that, given this discrete inner product, the 

representation is not continuous. As a consequence, the variables bg and b have no well- 
defined operator counterparts. Instead, one represents their complex exponentials, and 
that describe the holonomies of the connection: 


|y, — 


n ± 2, Ae ± 


2Afl 


e^^^\v,Xg) = \v±2,Xg). 


Note that, with these definitions. 


n}, Ag] = Ag}, 


( 2 . 2 ) 


(2.3) 


where ba denotes either bg or b. We also observe that the complex exponentials when 

acting on the above basis states, produce a shift on the quantum label v that is state- 
independent. This is the main reason motivating the above set of variables. Nevertheless, 
the translation on the anisotropy variable produced by is not constant. This is one of 
the features underlying the complications when solving the dynamics of the model. 

Besides, we take a standard Schrodinger representation for 0, with Hilbert space L^(M, d(j)), 
and momentum operator = —ihdfj,. Then, the Bianchi I term in the Hamiltonian con¬ 
straint is promoted to the following symmetric operator [ir 


fe2 

Gbi = Gfrw- Q —(G0 -|- 0H) ; Gfrw = 


37i^^^2 , § 
8 ^ 2 ■ 


Here, G is the Newton constant, Q = Qg — Q, with 

Q = \/i^ sign(n)sin(&) -|- sin(6)sign(n) \/|^, 


(2.4) 


(2.5) 


and Qg is dehned in a similar way, replacing sin(6) with sin(66i). These operators H and Qg 
represent, respectively, the classical functions 2nsin(6) and 2nsin(6e). The operator Gfrw 
is the Hamiltonian constraint of a flat FRW model coupled to a massless scalar 0. Then, 
the constraint operator Gbi of this LRS Bianchi I sector can be seen as an FRW term plus 
a contribution that accounts for the anisotropies. 

This Gbi operator has some remarkable properties thanks to the symmetric factor ordering 
chosen for H and Qg 13, It decouples the basis states |n, Xg) with n = 0 and/or A^ = 0 
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from their orthogonal complement, and it does not mix states with positive values of v 
and/or \q with states with negative values of those variables. This allows us to restrict our 
study to the subspace spanned by states |n, Xg) with, e.g., n, Xq G M'*'. It is now convenient to 
introduce a relabeling of the basis states in this subspace as |n. A) {y G M+, A G M), where 
A = ln(A6i). Moreover, this subspace further splits into separable superselection sectors 
under the action of Cbi- On the one hand, this action preserves all the subspaces spanned 
by states |n,A) with v belonging to the semilattice of step four Cf = {e + Ak,k G N} 
determined by the initial point e G (0,4], Notice that e is the (strictly positive) minimum 
allowed for the Bianchi I volume v in the considered sector. On the other hand, there are 
also superselection sectors in the anisotropy variable A. One can show that a state |u. A*) 
is related by the iterative action of the constraint only with states with A = A* + A^, where 
Ae belongs to the (countable and dense) set We dehned as 

I 2 ; In (- — ^ k^\n ; CeN, z e Z ii e > 2, z = ^\i e <2 \ . (2.6) 

I V ^ / m,nGN V ^ ^ / J 

Therefore, for the homogeneous gravitational part of the system, all our future analysis 
will be restricted to any of the sectors spanned by basis states |u,A) with v G and 
A = A* + Ae, with Ae G W^. 


B. Fock quantization of the inhomogeneous sector 


In the hybrid strategy, a Fock quantization is adopted for the non-zero modes of the 
inhomogeneous helds, expressing them in terms of a suitable set of annihilation and creation- 
like variables, that are represented as operators acting on the corresponding Fock space. 
Actually, the selection of a specihc quantization is possible thanks to the existence of a 
privileged choice of Fock representation for both the gravitational waves and the matter held, 
in the totally deparameterized Gowdy model. This representation is the unique one, up 
to unitary equivalence, that admits a unitary implementation of the dynamics and whose 
vacuum is invariant under rigid rotations in 0, which is the gauge symmetry of the reduced 
system (see [l5|). This class of unitarily equivalent Fock representations requires a particular 
choice of conhguration variables for the gravitational waves and for the matter held, that 
involves a homogeneous rescaling of those helds. The class contains the representation that 
would be most natural for massless free helds (that is, with annihilation and creation-like 
variables chosen as if the frequency of the modes did not include mass terms). 

With this result in mind, we adopt this “massless” Fock quantization for the non-zero 
modes of the inhomogeneous helds [ij. In this representation, we will denote the annihi¬ 
lation operator associated with the mode m E Z — {0} by dm\ where a = ^ denotes the 
gravitational held and a = (p denotes the matter held, both conveniently rescaled as stated 
above. The Fock space admits a basis of n-particle states, |n^,n‘^), where n“ denotes the 
inhnite collection of occupation numbers G N in each non-zero mode m of the held a. In 
this representation, the momentum constraint Ce that generates rigid rotations in 6 restricts 
the occupation numbers of the u-particle states so that T3 


meN+ 


m 




= 0. 


(2.7) 
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C. Hamiltonian constraint of the hybrid quantum Gowdy model 


The operator which represents the Hamiltonian constraint Cq in this hybrid quantization 


IS 


Co = C 




FRW 


(H0 + 0H) 


2TiGh^ 

~T~ 








( 2 . 8 ) 


Here, (3 = [Gh/( 167 r^ 7 ^A)]^C jg ^ constant that depends on some parameters of the loop 
quantization (7 is the Immirzi parameter 17j and A is the gap in the spectrum of area 
eigenvalues i), ^0 is the free-held energy of the non-zero modes, 


E (2.9) 

m€Z—{0} 

Hi is a self-interaction term, 

Hi= , ( 2 . 10 ) 

mEZ—{0} 

and D represents the product of the volume by its inverse, which is regularized in LQC j^. 
This product differs from the identity only in the region of small volumes: 

D\v) = D{v)\v), D{v) = V (^Vv -f 1 — ^/\v — l|j . (2.11) 


The last two terms in fl2.8p form the operator Cnh that represents the inhomogeneous term 
Cinh of the constraint. Given its action on the |n,A)-part of the states, the Hamiltonian 
constraint Cq preserves the superselection sectors introduced in Subsec. Ill A1 Actually, it is 
densely dehned on the kinematical Hilbert space dehned by the tensor product of L^(M, dip) 
with the completion of the linear span of |n,A,n^,n‘^) with respect to the inner product 


A', 


n'b A, n^, 




where we recall that v G Aj" and A = A* -|- A^, with A^ e We,. 


( 2 . 12 ) 


III. APPROXIMATING THE HAMILTONIAN CONSTRAINT 

The action of the Hamiltonian constraint Cq on the kinematical Hilbert space is quite 
complicated, a fact that makes very difficult (if not impossible) to analytically hnd the 
physical states of the system, those that solve the equation = oJl This is due to the 

presence of two terms: The anisotropy term containing ClQ + QQ and the interaction term 
containing e~'^^DC i^DH i. On the one hand, as far as the homogeneous sector is concerned, 
neither the anisotropy operator ClQ + 0f2 nor the operator D commute with the FRW 
operator Therefore they cannot be diagonalized simultaneously. This fact, together with 


^ Here |'I') generally stands for a generalized state. Alternatively, we can understand our constraint as an 
equation of the form ('I'ICq on generalized states on the dual of the domain of Cq. 













the u-dependent translations that the anisotropy operator produces in A [see fl2.2p ]. makes 
the resolution of the constraint in the homogeneous variables v and A extremely hard. On 
the other hand, concerning the inhomogeneous sector, the interaction operator Hi, which is 
densely defined on the linear span of the n-particle states, does not act diagonally on them 
(in contrast with Hq). Specifically, it creates and annihilates a pair of particles in every 
mode, notably complicating the resolution of the constraint for the inhomogeneous degrees 
of freedom. 

Owing to these obstacles, in what follows we will carry out approximations that allow 
us to disregard those problematic terms when acting on certain families of states. This will 
lead to a much simpler constraint operator, possessing some solutions that can be regarded 
as approximate solutions to the full Gowdy model, and which will be studied in Sec. HVl 

Let us start by reviewing the spectral properties of the operator (2^. It is essentially self- 
adjoint, with an absolutely continuous, non-degenerate, and positive spectrum [l^. In each 
of the superselection sectors, with support of n in a semilattice £+, the delta-normalized 
eigenstates of (2^ with eigenvalue G M’*', 


K) = 


vdCt 


»|f). 


(3.1) 


provide a resolution of the identity. The eigenfunctions can be chosen to be real and present 
a feature that will prove essential for our future approximations, namely, that when p S> 10, 
ep(v) is exponentially suppressed for v 5, p/2. On the other hand, for v 3> p/2, these 
eigenfunctions are oscillatory (see, e.g., [uj). The exponential suppression of the region 
n < p/2 is a characteristic of the quantum geometry effects in the context of LOG, and it 
is at the root of the occurrence of a quantum bounce in the loop quantization jig. In the 
proximity of this bounce, the LQG phenomena alter signihcantly the gravitational behavior 
with respect to the predictions of General Relativity, invalidating the expectations based 
on this latter theory [l^. Actually, this bounce resolves the cosmological singularities. 


Moreover, it persists in the presence of anisotropies [19| and of inhomogeneities [20 


Let us now consider states whose homogeneous gravitational contribution has the form 


'‘'■SJ. 


dcu g{u!. A* + Ag) jcu. A* -|- A^), 


(3.2) 


AeSWe 


' Spc(a;) 


with 


g{uj,A) = N{uj)f{uj,A), 


f{u. A) = e 2-^^ 


(3,3) 


Here, 


= In ( IH- 

'^m 


(3.4) 


with Vm £ a certain value of the volume such that Vm S> 10. Besides, as is a free 
parameter characterizing the width of the Gaussian /(cu. A) (together with g^), and Spc(a;) 
denotes the spectrum of some operator Cj defined on the homogeneous and isotropic geometry 
part of the kinematical Hilbert space (that is, the space spanned by the states |u)). We 
will assume Cj to be (essentially) self-adjoint. Let us notice that its spectrum might be 
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continuous, discrete, or even a mixture of both types. Nevertheless, we will formally denote 
the spectral resolution of the identity provided by a) as I = dci;|ci;)(a;|. States \Q) 

of this kind were studied in [Tl[, though only in two particular cases: First when cD is a 
constant operator, and then when uj = v. Moreover, keeping in mind the key properties of 


the states in [ll| needed for the different approximations introduced in those references, we 


will for the moment assume that the profile in the p-representation. 


9{P,^)= dw^(w,A)e^(p), (3.5) 

J Spc(i.j) 

is highly suppressed in the region with p ^ Pm = ‘2vm- Here we are denoting by e^(p) 
the wave function of the state ja;) in the p-representation, e^(p) = {e^p\oj). Owing to the 
exponential suppression of e^(n) in the region v < p/2 discussed in the previous paragraph, 
a prohle g{p,A) highly suppressed for p < p^ in turn implies that the corresponding prohle 
in the ^-representation, g{v,A) = dpg(p,A)ep(v), is highly suppressed for v 5 , Vm- 

In the following, we will analyze under which conditions the approximations of [1I| extend 
to these generalized states with homogeneous prohles characterized by (I3.3j) . Later on, in 
Sec. m we will explain how to construct approximate solutions with the desired properties. 


A. Approximating the anisotropy term 


In order to deal with the anisotropy operator ClQ + it is essential to notice that, 
when acting on the considered states \Q), and owing to the suppression in the region v <Vm 
that we are assuming, the only contributing v-dependent shifts on the anisotropy variable A 
that its action produces are not (significantly) bigger than the scale ll|. Recalling that 


the considered profiles g{u, A) are Gaussian-like in the anisotropy variable with width given 
by qe/cTs, if that scale is much smaller than their width, namely if as -C 1, then they can be 
extended (from the anisotropy superselection sector) to a smooth function in A so that, for 
contributing shifts Aq < q^, 


g(ui, A + Ao) ~ g(w, A) + Ao9a9(w, A), 


(3.6) 


for all UJ in the support of g(uj, A). If one considers the ^-representation of \Q), this in turn 
implies that g{v, A + Aq) — g{v, A) -|- A^df^giy, A) for all v. 

For such states, the operator G© -|- ©G approximately factorizes as —2VlQ'. Here, G is 
dehned like the geometry operator G in fl2.5p . except for the substitution of the conjugate 
pair of variables (w, 5) by (n/2,2&). On the other hand, ©' is the discretization of the first 
derivative —did^ at the scale q^: 


0'|A) = i— (|A + qs) - |A - qs)) . 

Qe 


(3.7) 


For further details. 


The proof of this approximation follows exactly the same steps as in [11 
we refer the reader to that work. 

Note that q^ G We so that G©' preserves the superselection sectors where ClQ + QQ is 
defined. Furthermore, the action of 0' leaves invariant the lattices of constant step q^ of 
the form = {A' -|- nge; u G Z} contained in those sectors, with A' — A* G Wg. Using 
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this result, we can further restrict the study of the anisotropy variable just to the subspace 
spanned by states |A) with support in any of those lattices. 

We will now show that the approximate anisotropy operator given by —2QQ' can be 
disregarded when acting on the considered states \Q) in comparison with First of all, it 

is easy to check that our definition of the operator hi is equivalent to the following action on 
the basis states: 


fi|n) = i [y-{v)\v - 4) - y+{v)\v + 4)], (3.8) 

with 

S±(«) = 1 + ± (3,9) 

Employing this, we get that the approximate action of the negative of 110 + 011 on our 
states is given by 

(n,A|2il0'|^) = — {y+{v) [^(n + 4,A + g0 -^(n + 4,A-g0] 

Qe 

-y-{v) {g{v - A, A + q,) - g{v - A, A - q^)] }. (3.10) 

Now, given the form fl3.3p of the profiles, a straightforward calculation leads to 
2 

— [giv ± 4, A + - g{v ± 4, A - q,)] 

Qe 

= ——e 2 " f dci; 5 f(a;, A)e^(n ± 4) sinh [ — [A — A(a;)] ] . (3.11) 
*?£ J Spc(L»;) V / 

Here, e^(n) = (nlcn) denotes the wave function of the state Icn) in the n-representation. The 
function g{u,A) only contributes when crs[A — A{u)]/qe is 0{1). Therefore, we get 

^ 4 miT n (7 

— [g{v ± 4, A + - g{v ± 4, A - ~- -e~^g{v ± 4, A) x 0(1). (3.12) 

Qe Qe 

Hence, recalling (I3.10p . we have just obtained that the action of the Bianchi I anisotropy 
operator on these states can be approximated as follows, when one stays in the sector 
V ^ Vm ^ W under consideration: 

{v,A\Cie + 00|^) ~ (n, Alfil^), (3.13) 

Qe 


up to a factor 0{1) on the right-hand side. Taking into account that the action of 11 on \Q) 
is of the same order as that of H, since these operators are completely analogous except for 
the magnitude of the shifts that they produce in n, we can expect the right-hand side of 
fl3.13p to be negligible compared to (n,A|ll^|^) if 8 | sinhcTsI exp (— cr^/2) qeVm- Here, we 
have used the fact that, in the region with v > Vm^ 10 , 


0^0) = —y^+{v)\v -f 4) -|- 2v^\v) — y _(n)|n — 4), 


(3.14) 
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with 


y++{v) = y+{v)y+iv + 2), y—{v) = y++{v - 4), (3.15) 

1 + sign(u + 2) / 

y+{'v) = -^-^^^(^ + 2 ), ( 3 . 16 ) 

and so the coefficients of {v, A|r2^|^) are 0{v‘^). Equivalently, since for Vm ^ 10, ~ 2/vm, 

we can rewrite the above condition as 

4|sinh(Js|e ^ <C 1. (3-17) 

Recalling that we have already required that cTs <C 1 in order to treat g{v,A) as smooth 
in A and approximate it by its Taylor expansion truncated at the first derivative, the new 
condition fl3.17p turns out to be trivially satished in the considered case. 

This concludes the proof that the action of the anisotropy term on states with profile of 
the form fl3.3p can be disregarded in the Hamiltonian constraint provided that the condition 
(Ts 1 is satished and that g{p,A) is highly suppressed for p < pm- This was actually 
expected, since the Gaussian-like prohles that we are considering for the anisotropies are 
reasonably centered on trajectories with vanishing momenta 0'. 


B. Approximating the interaction term 

Let us analyze now the interaction term e~‘^^DCl^DHi. Before discussing whether this 
operator can be neglected in the Hamiltonian constraint when considering states \Q) with 
prohles of the form fl3.3l) . it is worth noticing that, owing to their suppression in the sector 
v ^ Vm with Vm ^ 10, scctor on which D just acts like the identity according to 02.111) . we 
can make the approximation; 




(3.18) 


Thus, provided that the content of inhomogeneities of our states is reasonable, we can 
disregard this interaction term if the action of on them is negligible compared to the 

rest of contributions in the Hamiltonian constraint. This action is given by 

Note that the hrst exponential of the right-hand side is bounded from above by the unit. 
Therefore we see that, if we choose the peak of the Gaussian-like profiles such that A(a;) is 
much bigger than both 1 and ql/ul for all u in the support of N{uj), then this contribution 
will be negligible compared to In conclusion, we can disregard the interaction term 

in the Hamiltonian constraint of the Gowdy model for states with the considered profiles 
if it is satisfied that A(a;) S> max (1, g^/a^) for all values of co in the support of the state 
(where max stands for the maximum). 


ie>= E 




'Spc(Lj) 


dcu N{u)e 


2<jf 


A-A(a;)+ 


2qt 


J e 


-2A(i.j)+ 


21? 


■0V)- 


(3.19) 
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C. Approximating the free term 


There is one more approximation that can be made in the Hamiltonian constraint Cq 
on states \Q) with the prohle (13.3p . This is an approximation for the term containing the 
free contribution of the inhomogeneities, namely (up to constants) The idea is to 

restrict the Gaussian-like prohles in g{uj, A) to be sharply peaked at A(a;), that is, to require 
its width qe/us to be much smaller than the unity. If this is so, then 

~ (3.20) 


Indeed, the support of a Gaussian is approximately its width, and in our case this support 
corresponds to A such that |A — A(a;)| < qe/ag- Therefore, the non-negligible contributions 
of the action of on these states will take the values 


2A 2k(u))+2a^ 

e = e ^ 


a e [-1,1]. 


(3.21) 


Thus, if qe/c^s -C 1, this value will be essentially equal to that taken by the operator 
on our state (term by term in its spectral decomposition). Then, this approximation will 
be consistent with disregarding the anisotropy term if and only if the parameters that 
characterize the width of the Gaussian-like prohle are such that -C 1. Let us 

conclude noticing that, in this way, by considering states with prohles for the anisotropy 
sharply peaked at a function of some operator of the homogeneous and isotropic geometry, 
one ends up mimicking a contribution in the approximated Hamiltonian constraint that is 
given precisely by that very same operator, in the sense of fl3.2Up . 


IV. APPROXIMATE SOLUTIONS TO THE GOWDY CONSTRAINT: 
MODELING MODIFIED FRW COSMOLOGIES 


In the previous section, we have shown how the full Gowdy Hamiltonian constraint Cq 
can be approximated by the operator 


r' -r , _ 

^aoD “ LfRW H - 75-^ -^0 — 


8 2 (3 ° 


app ^rrxvv I p 

when acting on quantum states whose homogeneous gravitational part is given by 

da;g(a;, A)|a;, A), 




(4.1) 


(4.2) 




Gpc(a;) 


with g{u,A) dehned in (13.3p . provided that the following conditions are satished: 

i) In the p-representation, g{p, A) has to be highly suppressed for p < pm, with pm S> 10. 
This in particular implies that g(n. A) is exponentially suppresed for v < Vm = Pm/‘2. 

ii) A(ci;) 3> 1 for all u in the support of N{u). 

hi) qe <^as 1. 

In what follows, we will consider solutions of this approximate constraint that are in turn 
approximate solutions of the full Gowdy model. Later on, a possible physical interpretation 
of those solutions will be given. 
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A. Construction of solutions 


Let us search for solutions of the approximate constraint Gappl^) = 0 by considering 


states 


1^)=/ [ da; ^(P0,w,A,n«,n^)|p0,a;,A,n«,n'^) (4.3) 

J-oc JSpcM 

A' ’ 


with wave function of the form 


'L(p0,a;,A,n^n‘^) = /(a;, A)A(a;,p^, n^), (4.4) 

where /(a;, A) is given in fl3.3l) and satishes by construction our condition iii) above. The 
sets of occupation numbers and determine the eigenvalue of Hq (which acts diagonally 
on the n-particles states): 

Ho = Ho{n^,n‘^)= ^ \m\{nl^ + n‘^). (4.5) 

mGZ—{0} 


Here, we have generically included in N{u)) the dependence of the wave function on and 
on these occupation numbers. This function N{u,p^, n^, n"^) should be chosen in such a way 
that the content of inhomogeneities is small, so that the approximation of disregarding Ci 
holds, and such that the momentum constraint (12.Th is satished. 

Since p^ and Hq are Dirac observables of this approximate constraint, and therefore p^ 
and Hq are constants of motion, solving the approximate constraint C^pplT) = 0 on the 
above states is equivalent to solve (TapplT) = 0 in each eigenspace of the two considered 
Dirac observables, with (Tapp being the operator 


Capp — ~ 




2 


27r(Th2 


f3 


(4.6) 


Note that this constraint operator only acts on the homogeneous and isotropic part of the 
Hilbert space, namely that of the FRW model. However, let us emphasize the fact that, to 
each state T, there corresponds a (collection of operators) (Tapp, as this operator depends 
(besides as on p^ and Hq) on 

Keeping in mind that we are interested in states that are in turn approximate solutions 
of the full Gowdy model, the wave function 04.41) must be such that these states satisfy 
conditions i) and ii) above as well. In order to fulhll these conditions, and motivated by 
the strategy employed in Tl[, we will restrict all considerations to states with a peak of the 
Gaussian A(a;) such that the resulting operator A(a;) in 04. 6 p is defined in the following way: 


Ao, if n < no, 

A(a)o), if n > no. 


(4.7) 


for certain no > Vm- Here, Ao is a constant and Do is the restriction (via projection) of D to 
the sector of the homogeneous and isotropic part of the Hilbert space with n > no. Explicitly, 
Do = PDP, where P is the projector on the linear span of |n) with n G = Cf D {n > no}. 
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Note that Vq does not necessarily belong to the semilattice considered. For future 
reference, let us denote by Vi the lowest end of so that Uq < + 4. Furthermore, 

we will focus on operators Cj such that 

g2A(^o) = g2ft(«) ^ 

where the function h{v) varies sufficiently smoothly so as to allow us to make the approxi¬ 
mation h{v ± 4) ~ h{v) in the region v > uq- Besides, we assume Op to be a positive (and 
thus self-adjoint) operator dehned on the linear span of |u) with v G with a quasi-local 
action on this basis of the generic form 

K 

Op\v) = ^[f^{v)\v + 4:k) + f^{v)\v-4:k)], K < oo. (4.9) 

k=0 


Here, /o’ = /o = /o and (v) = 0 ii v — Ak < vq, so that Op is indeed dehned on the sector 
V > Vq- The case with Op = 0 was the one studied in IfJ. For conciseness, we now assume 
that /^(u) is (strictly) positive for v > Vq. As shown in Appendix |3, states T that solve 
the constraint equation C'appl^') = 0 can be determined from the equation 


N{v,p^,n^,n‘^)C^pp\v) = 0, (4.10) 

V&C+ 


with 


N{v,P4„n^,n‘^)= [ dw A^(a;,p^, n^, n‘^)e^(u). (4.11) 

J Spc(a;) 

We will now analyze how we can construct solutions to fl4.10l) with the desired properties. 
First of all, condition ii) is automatically satished for these states if Aq S> 1 and h{y) S> 1 
for all u G Concerning condition i), for v < vq the operator A(a)) is simply a constant 
and the constraint equation reduces to an eigenvalue equation for the FRW operator (1^, 
given by H 

(4.12) 

Therefore, for all v < vi solutions are of the form 


where we recall that is the wave function in the u-representation of the eigenfunction of 
with eigenvalue p^. Hence, condition i) for these states is satished, for instance, if one 
restricts the function ,vO) to have support on the region with p^ > p™, where we 

choose ^ \/757rGh. That this is true follows from the positivity of the last term in 
fl4.12p . Indeed, in that case, solutions will have signihcant contributions only for 


P > 


:P4> > 


VSirGh vSttGK 


p™ = » 10. 


(4.14) 
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The identity in this formula is just a dehnition of the scale pm used in our construction of 
states. Let us note that this scale (and therefore the corresponding value of Vm) has been 
dehned in an intrinsic way, in terms of the conserved momentum of the homogeneous scalar 
held. 

In the region with v G equation (I4.10j) leads to the difference equation [dropping the 
dependence of iV(n,p0, on and the occupation numbers in order to simplify the 

notation] 
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SuGh? 3/9 

K 


N{v) + y++{v)N{v + 4) + y _ {v)N{v — 4) 


[/^(^“4/c)iV(n-4/c)x£p(n-4fc) + 4 {v + Ak)N{v + Ak)] =0. (4.15) 


fc=i 


Here, 


Xc^iv) 


1, if n e 

0, ifn^£f. 


and we have dehned the function 

A(n) = ^ln [e2M-)+2/o(n)]. 


(4.16) 


(4.17) 


The diherence equation fl4.15l) involves coefficients N{v) evaluated on K points above the 
considered one, n, in the semilattice. In particular, this happens at the matching point with 
the solution of constant Aq, i.e., when v equals the lower end vi of Owing to this fact 
and that, at this matching point, we know data only for values of v smaller than or equal to 
it, it is not difficult to realize that, in order to be able to hnd an approximate solution to our 
equation without ambiguities, we must impose certain requirements on the operator Op. For 
instance, to arrive at acceptably smooth solutions, we can require that the functions f^{v) 
with A; 7^ 0 be negligible in a neighbourhood above vq, namely, at least at all points in the 
interval X = [uq, "^0 + 8K — 4). Indeed, if this is so and after disregarding the contribution of 
those functions, then fl4.15l) approximately gives us in a deterministic way the K — 1 values 
of the function N(v) from n = ui + 4 up to n = Ui + 4K — 4, when supplemented with the 
input data 


ep(p„Ao,Ho)(^i) and (4.18) 

coming from the imposition of the constraint in the region v < Vq < Vi. We note, in 
particular, that to fix A^(ni + 4K — 4) in this manner we need to ignore the functions f^{v) 
at least up to Vi + 8K — 8, for any possible value of vi G [uq, "^0 + 4), as we have assumed above 
that it is indeed the case. The K coefficients computed in from N{vi) to N{vi + 4K — 4), 
can then serve as initial data in order to uniquely fix the rest of coefficients. The procedure 
to do so is to consider again the full constraint equation fl4.15p evaluated at points v > Vi 
without neglecting the functions f^{v). Introducing the value of those K coefficients, the 
constraint for v = Vi completely determines the next coefficient N{vi + 4K), and so on and 
so forth for v > Vi + 4. With this method, the whole solution can be constructed, at least 
approximately. 













16 


Let us notice that, if the neighbourhood X contains n > 2K — 1 points of then the 
number of coefhcients that can be obtained by ignoring the functions f^{v) in fl4.15l) is 
n — K > K — 1. In this situation, the approximate solution that has been obtained can be 
improved by iteration at the n — 2K + 1 points that are just above Vi. This can be done by 
considering again the constraint equation (14.151) evaluated at the point just below the one 
that we want to improve, taking into consideration the corrections given by the contribution 
of the approximate coefficients of the K — 1 larger nearest points, and (up to) the K + 1 
smaller ones in X^, all of them multiplied by the corresponding functions f^{v). 

Finally, it seems natural to impose continuity of the approximate constraint when v = vi. 
From the given construction of the solutions, it is straightforward to check that, for this 
continuity to be obtained, it suffices to £x the constant Aq = A(ui), with A(ui) ::§> 1. 

We refer the reader to Appendix [B] for additional details on the suppression of contribu¬ 
tions with p < Pm in the solutions that we have constructed. 


B. Perfect fluids and geometrical corrections 


In the previous section we found some approximate solutions to the hybrid Gowdy model 
that, besides, effectively obey a dynamics dictated by the constraint 


i-^app Gppfw 


27rGh2 r 

~T~ 


,2h{v) ^ Q 


Hn 


(4.19) 


for sufficiently large volumes {v > vq). We will now discuss how this constraint can be 
understood as the one corresponding to an isotropic flat FRW model coupled to different 
types of perfect fluids coming from the term and geometrically corrected by the term 

Op, which can be interpreted as arising from homogeneous curvature-like terms or higher- 
derivative contributions in the gravitational action. Indeed, as discussed in lH, if we choose 


h{v) =lnQ5^ 


(l-W)) 


11/2 


(4.20) 


the dynamics of the constructed states mimics that typical of a content of different perfect 
fluids with equations of state given by p = we, where p and e denote respectively the pressure 
and the energy density of the corresponding fluid. Here, w runs over as many parameters as 
wanted (one for each different perfect fluid), provided that ta < 1 so that the approximations 
done in Sec. m are valid. This upper bound for w allows for physically interesting couplings 
such as dust, radiation, and a cosmological constant, obtained by setting w = 0, w = 1/3, 
and w = —1, respectively, in our formulas. Note that the conditions on the peak of the 
Gaussian-like prohles that are needed for our approximations to hold are then automatically 
satished if Vq S> In addition, we can also treat the case tc = 1 as corresponding 

to a massless scalar held contribution, coming from the exponential term in the dehnition 
fl4.8p . that can be included in the homogeneous held (j) of the FRW constraint [Tl| (just by 
redehning the latter). 

Regarding the new effective term Op, its action on the kinematical basis provided by 
|u), given in fl4.9p . can be constructed, for instance, from sums of powers of possibly 
multiplied by smooth functions of v. Hence, taking into account that the FRW operator 

fully characterizes the curvature scalar R of a flat FRW universe [^, the contribution 
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of Op in this modified constraint Capp can be seen as a term corresponding to additional 
curvature-like terms correcting R in the gravitational action. Alternatively, it is also possible 
to interpret it as discretized higher-derivative terms. These kinds of terms are some of those 
that one would expect to appear in certain /(i7)-theories and other modihed theories of 


gravity (see, e.g., jl2l]h 


To sum up, we have seen how some approximate solutions of the Gowdy model, that is 
genuinely anisotropic and inhomogeneous, can effectively behave as solutions (also approx¬ 
imate in general) of the Hamiltonian constraint of a flat FRW model coupled to different 
types of perfect fluids and with geometrical corrections similar to those of modihed gravity. 
It is worth clarifying that, despite of the dynamical behavior proven for these states with 
respect to the constraint of the system (namely, a homogeneous and isotropic effective con¬ 
straint), their Gaussian-like prohles are not peaked on isotropic trajectories of the classical 
model, but generically on trajectories that are very anisotropic. Glassically, isotropy implies 
the relation 3A = In (n/2), that ensures that Xq = coincides with the geometrical average 
scale factor of the model, given by the cube root of the volume, up to proportional¬ 
ity factors 1^. The Gaussian-like wave functions fl3.3p . however, may be peaked on many 
possible trajectories, determined by quite general functions A(a)) of a variety of homoge¬ 
neous and isotropic operators u. As we have pointed out, in general, these trajectories do 
not correspond to isotropic and homogeneous solutions of the classical Gowdy system. At 
the end of the day, it is in the collective behavior of the anisotropies and inhomogeneities, 
together with the quantum effects of the loop quantization of the geometry, where one hnds 
the ultimate reasons explaining the approximate dynamics of FRW-type, with curvature- 
like and higher-derivative corrections, that the considered states display. In this sense, the 
geometrical modihcations to the FRW dynamics obtained for these states may be regarded 
as arising from the underlying quantum theory (both from the loop representation and from 
the characteristics of the states). Finally, it is worth remarking that such an effective de¬ 
scription starts to apply only when one reaches the volume vq- For smaller values of v we 
get instead an effective dynamics which is just that of an FRW model coupled to a ho¬ 
mogeneous massless scalar held. Nonetheless, once the epoch with geometrically modihed 
dynamics and perfect huids is reached, that regime holds indehnitely, for all n > Uo) by the 
very construction of our states. 


V. CONCLUSIONS 

We have investigated the construction of approximate solutions in the hybrid quantum 
Gowdy model with three-torus topology, linear polarization, LRS, and a minimally coupled 
massless scalar held [l^. More specihcally, we have managed to construct approximate 
physical states of this inhomogeneous model that in turn are also (approximate) solutions to 
the Hamiltonian constraint of a homogeneous and isotropic hat FRW model with corrections 
that can be interpreted as curvature-like or higher-derivative terms. The present work 
signihcantly generalizes the results of 0 , where we already provided approximate solutions 
of the Gowdy cosmology that resemble (as far as the Hamiltonian constraint is involved) 
those of a hat FRW universe with a massless scalar and a perfect huid. Based on the 
approximations developed in those previous papers, and extending the analysis carried out 
there, now we have constructed approximate solutions that behave as those of a geometrically 
modihed hat FRW containing diherent types of barotropic perfect huids with equation of 
state characterized by parameters tc < 1. 
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Our results show how some specihc quantum solutions of inhomogeneous models, in this 
case the Gowdy model, can behave dynamically as solutions of flat homogeneous and 
isotropic cosmologies with a particular kind of homogeneous and isotropic matter content, 
and even with homogeneous and isotropic geometrical modihcations that can be regarded 
as higher powers of the curvature or higher derivatives. It is worth emphasizing that those 
solutions are far from being genuinely isotropic and homogeneous. Their anisotropies and 
inhomogeneities are not negligible, as one could show, in principle, by measuring on those 
states generic quantum observables beyond homogeneity and isotropy. Despite of that, 
it is remarkable that those states behave in such a way that they lead to effective terms 
in the Hamiltonian constraint which are characteristic of a flat FRW model, and more 
specihcally an FRW universe in the presence of perfect fluids and geometrical corrections. 
In particular, the dependence of the peak of the considered Gaussian-like prohles on a 
homogeneous and isotropic geometric operator (of a certain type, though quite general) 
turns into the appearance of that operator in the effective constraint. This phenomenon, 
that strongly depends on the specihc choice of the considered family of states, emphasizes the 
fact that effective descriptions generally depend on the particular set of states under analysis. 
This is an idea which is attracting increasing attention lately (see, e.g., [1^ for discussions 
on other types of inhomogeneous quantum states for which simple homogeneous descriptions 
are obtained). In fact, the effective dynamics attained here can be understood to arise from 
the quantum correlations existing in the considered states between the different sectors of 
the homogeneous Hilbert space, namely the set of states studied here presents prohles with 
a specihc mixed dependence on the variables of the homogeneous phase space. Besides 
the mentioned correlations, two key properties of the constructed solutions lay behind this 
interesting behavior: They display a negligible momentum of the variable that measures 
the anisotropy and they experience a negligible coupling between the homogeneous sector 
and the self-interaction of the inhomogeneities. It is because of these properties that one 
can disregard the most problematic terms in the Hamiltonian constraint, arriving then to 
a much simpler constraint operator corresponding to a modihed hat FRW model coupled 
to perfect huids. It is worth mentioning that this approximation is consistent thanks to 
the quantum geometry ehects introduced by LQC, that in particular are responsible of the 
exponential suppression of the eigenstates of the FRW geometry at small volumes. These 
quantum ehects, together with the collective behavior of anisotropies and inhomogeneities, 
produce departures from the typical classical behavior predicted by General Relativity in 
the high-curvature regime, namely, around the cosmological singularity, where the role of the 
anisotropies would have been very relevant. This is the reason why the ehective dynamics 
obtained for these families of states diher in those regions from the classical solutions of the 
model. 

Our analysis sheds light as well on the generality of our approximations and on the extent 
to which they are robust. We may expect that these approximations continue to hold as far 
as one considers states with negligible contributions of the sector p < pm (with pm 3> 10) 
for the FRW geometry operator, and anisotropies peaked around a certain trajectory on the 
homogeneous phase space, with a large value of the peak, such that two more conditions are 
satished. First, the momentum of the anisotropy variable must be small, with a dispersion 
negligible compared to pm, that when squared can be considered a lower bound on the 
FRW geometry operator acting on our states. Second, the dispersion in the anisotropy 
variable (dual to that in the momentum of such variable) must be much smaller than the 
unity. As we explained, the hrst condition allows us to neglect the anisotropy term, while 
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the second, together with the requirement of a high value of the peak in the anisotropy, 
permits to neglect the term that contains the self-interaction of the inhomogeneities, and 
approximate the term with their free contribution using FRW operators. These arguments 
indicate that Gaussanity on the anisotropies is not strictly necessary. They also tell us that 
our approximations may be considered stable within a certain sector of quantum states, as 
long as we do not compromise any of the required properties. 

Another issue that may be helpful in understanding the bases of our approximations is 
the existence of asymptotic limits in which they arise naturally. That such limits can be 
found whithin certain schemes follows from the requirements imposed on our states at the 
beginning of Sec. HVl Assuming states peaked on large values of the anisotropy variable, and 
recalling that, for large pm = ‘^Vm, one has — ^1 Pm-, we can proceed in the following way. 
First, we focus our attention on functional relations of the form = q^ with o<C< 1. For 
instance, we can take C = 1/2. Then, the necessary conditions for our approximations are 
all reached, with improving accuracy, in the asymptotic limit in which p^ tends to inhnity, 
in the sector of inhnitely large eigenvalues of the FRW geometric operator. 

Finally, it may be tempting to extrapolate the lessons learned here about the peculiar¬ 
ities of the effective dynamics associated with certain sets of states in order to build new 
avenues for the resolution of some of the open questions of the standard cosmological model 
with inflation. For instance, it may be worthy to investigate whether some elements of its 
phenomenology, such as the existence of a small but non-vanishing cosmological constant, 
or the origins of inflation, can be understood in terms of an effective description that arises 
from a global behavior in the quantum realm, accounting for a multitude of extra degrees 
of freedom, possibly inhomogeneous and anisotropic in nature. Besides, from a perspective 
which is beyond that of effective held theories, the results obtained in this work suggest 
the possibility that some of the corrections to the Einsteinian theory that are nowadays 
widely investigated in the context of the so-called modihed theories of gravity, may actu¬ 
ally be rooted in more fundamental quantum geometry effects, occurring in certain types of 
quantum states. 
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Appendix A: Validity of the constraint eqnation 

For each sector of constant eigenvalue of the Dirac observable and occupation numbers 
given by and (and hence also with constant eigenvalue of the Dirac observable 77o); 
and for all A G £^/, the constraint equation corresponding to Gapp can be written on the 
states fl4.4l) as 


' Spc{Lj) 


dw V(a;,P0,nAn‘^)Gappe 


A-A(il)) 


|u;) = 0. 


(Al) 
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Let us consider the operator H 







(A-A(w))^ 


(A2) 


Then, we can approximate the constraint by 

2 r — 12 

g ^ iV(n,p<^,n^,n‘^)C'app|n) = 0 (A3) 

^ 6 £+ 

if we can neglect the contribution of flA2p as a correction to Capp in this equation. If this is 
so, the profiles 


N{v,p^,n^,n‘^) = [ duN{uj,p^,n^,n'^)el{v) (A4) 

J Spc(tj) 

that satisfy fl4.10p provide indeed approximate solutions to our constraint. 

Let us hrst recall that, given two operators A and B whose commutator is negligible in 
comparison with one of them, the Baker-Campbell-Hausdorff formula implies (see, e.g., 0 ) 

e*Be-* = B+[A,B] + AMA B]] + ■ ■ ■ ~ B+[A,B], (A5) 


where we have ignored higher factor-order contributions, as the commutator [A, B] is as¬ 
sumed to be negligible. Applying this equation, we can rewrite our condition on flA2p as the 
same requirement on 



(A-A(ch))^(7app 


~ — (A 

(fe 


A(w)) (Tapp, A(w) 


(A6) 


On the right-hand side of flA6p . we have ignored again higher factor-order contributions. 
Taking into account that the only operator in (Tapp that does not commute with A(a)) is 
and that, regarding the behavior in A, the Gaussian-like prohles of our states only contribute 
when (Ts[A — A(a;)]/g£ is C>(1), the condition on the right-hand side of flA6p (ignoring again 
factor ordering and irrelevant factors) amounts to demand that 


^g- 2 A(*) 

Qe 



(A7) 


can be neglected if considered as a correction to the constraint (Tapp- 

Let us check that this condition is satisfied, and therefore that flA3p is approximately valid 
on the studied solutions. Since, according to fl4.7p . A(d;) is constant in the region v < Uoi 
the corresponding contribution to the commutator appearing in flA7p vanishes. Therefore, 
it suffices to analyze the similar commutator obtained by replacing u with cuq in the region 
V > Vo- Recalling definition fl4.8p for ujq, there are two terms to discuss in the commutator: 

and Op]. As for the first of them, the function h{v) has been chosen to satisfy 


Here, we replace square brackets with parentheses in the expressions involving A — A to avoid confusions 


with commutators. 



















21 


h{y ± 4) ~ h{v) for v > vq. Therefore, since the action of given by (I3.14p . essentially 
shifts the volume of the state in four units, up and down, it turns out that we can discard 
the anayzed commutator in our approximations. The second term requires a more careful 
analysis. On the one hand, the action of this commutator on |n)-states with v > V\ + AK 
yields 

K 

.odi-> = E 

k=0 

+ [y++{v)fk + 4) - y++{v - 4/c)/“(n)] \v - A{k - 1)) 

+ - 4) - y__{v + 4/c)/+(n)] \v + A{k - 1)) 

+ [y—{v)fkiv -4:) -y—{v-Ak)f-{v)] \v-A{k + l)) 

+ 16{kv + 2P)f^{v)\v + Ak) — 16{kv — 2P)f^{v)\v — 4/c)|. (A 8 ) 

These terms [and a fortiori the corresponding ones in flATIl ] will give negligible corrections 
to (Tapp on our states if Op is such that 

a) The functions f^{v) are smooth enough as to satisfy f^{v ± 4) ~ f^{v) in the con¬ 
sidered region v > Vq. 

b) The integer K is small enough so that y++{y ± Ak) ~ |/++(n), y _ {v ± Ak) ~ y _ (n), 

16kv <C and -C for all k < K in the sector with v > Vq- Using the fact 
that, according to their dehnitions fl3.15p and fl3.16p . the functions 2 /±±(n) are O(v^) 
in the studied sector, we can say that an integer K would be sufficiently small in this 
sense if it satishes K uo/lO. 

On the other hand, if one considers the action of the commutator [O^, 0f\ on the remaining 
sector of |n)-states, with v < vi + AK, other contributions appear different from those in 
flASp . This peculiarity occurs because is dehned on the whole semilattice Cf, whereas Op 
is dehned only on the restriction C^. These additional contributions are terms of the form 

K K 

'^y—{vi)f-{vi + Ak)\vi-A), '^y^+ivi - A)f+{vi)\vi + Ak). (A9) 

k=0 k=0 

Here, the hrst term accounts for the action of the commutator on all the states |ni -|- Ak) 
with 0 < k < K, and the second one includes the contribution of the action on |ni — 4) 
(recall that vi is the lowest point in £§). As a consequence, if conditions a) and b) hold, 
and 

c) the functions fki^i + Ak) and fki^i) are much smaller than the unit for 0 < k < K, 

then all terms under discussion will give negligible corrections to the constraint equation 
corresponding to (Tapp. Let us comment that this last condition c), for fc 7 ^ 0, was required 
as well at the end of Subsec IIV Al in order to be able to determine the solutions to the 
approximate constraint. Taking that into account, here we are just including a similar 
requirement on fo{vi). 




[y++{v)f+{v + 4) - y++{v + 4fc)/+(n)] \v + A{k + 1)) 
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Appendix B: Suppresion of the solution when p < pm 


In Subsec. II V A I we discussed the construction of approximate solutions to the constraint 
Capp. In particular, our analysis contained a region v < vq where the peak of the Gaussian 
in the anisotropies is constant. In that region, we made sure that the solution, determined 
by (I4.13p . possesses only non-negligible contributions for p > Pm by choosing pm as in (14.141) . 
One may wonder whether this statement is still true for the whole solution, built for v > Vq 
with Gaussian-like peaks that are not constant anymore. Although this is not directly 
granted, because it involves the spectral decomposition of the whole solution in terms of 
eigenfunctions of the operator we show here that the change from fl4.13l) to the new 
solution in V E (i.e., for v > vo) respects that all relevant contributions have p > pm, as 
required. 

Up to a global numerical factor, the constraint operator Gapp, given in fl4.6p . can be 
rewritten in the form F — where 


4p2 


16 


r = — 

SnGh? 3/5 


(Bl) 


For our discussion, the important fact that complicates the analysis in the region v > vq 
is that the FRW operator does not commute with F. Nonethelesss, when acting on 
solutions to the constraint, the relation = F holds. Besides, for the quantum states 
that were considered in Subsec. IIV A1 characterized by prohles with support on values of 
the scalar held momentum > p™, the action of the operator F is always greater than 

(multiplication by) p^ [see fl4.14p ]. i.e., (F — p^) is strictly positive, because so is the last 
term in flBip . Hence, on the sector that contains our states, the action of will lead only 
to contributions with p > pm, as we wanted to show. In fact, one can demonstrate that, on 
the considered sector, it is a good approximation to neglect the commutator of and F, so 
that it is acceptable to work assuming that they can be simultaneously diagonalized, a fact 
that supports the conclusion presented above. 
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